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We study how well magnetic models can be implemented with ultracold bosonic atoms of two 
different hyperfine states in an optical superlattice. The system is captured by a two-species Bose- 
Hubbard model, but realizes in a certain parameter regime actually the physics of a spin- 1/2 Heisen- 
berg magnet, describing the second order hopping processes. Tuning of the superlattice allows for 
controlling the effect of fast first order processes versus the slower second order ones. Using the 
density-matrix renormalization-group method, we provide the evolution of typical experimentally 
available observables. The validity of the description via the Heisenberg model, depending on the 
parameters of the Hubbard model, is studied numerically and analytically. The analysis is also mo- 
tivated by recent experiments [S. Foiling et al, Nature 448, 1029 (2007); S. Trotzky et al, Sience 
319, 295 (2008)] where coherent two-particle dynamics with ultracold bosonic atoms in isolated 
double wells were realized. We provide theoretical background for the next step, the observation 
of coherent many-particle dynamics after coupling the double wells. Contrary to the case of iso- 
lated double wells, relaxation of local observables can be observed. The tunability between the 
Bose-Hubbard model and the Heisenberg model in this setup could be used to study experimentally 
the differences in equilibration processes for nonintegrable and Bethe ansatz integrable models. We 
show that the relaxation in the Heisenberg model is connected to a phase averaging effect, which is 
in contrast to the typical scattering driven thermalization in nonintegrable models. We discuss the 
preparation of magnetic groundstates by adiabatic tuning of the superlattice parameters. 

PACS numbers: 37.10.Jk, 75.10.Jm, 05.70.Ln, 02.30.Ik, 



I. INTRODUCTION 

One of the most exciting recent events in physics 
has been the increasing overlap between two previously 
disjoint fields, quantum optics and condensed matter 
physics. This has become possible due to the enormous 
progress in cooling dilute bosonic and also fcrmionic gases 
down to temperatures where respectively Bose-Einstein 
condensation and Fermi degeneracy (temperatures well 
below the Fermi energy) are reached. 

A very attractive feature of this new class of exper- 
iments is that they provide the arguably cleanest real- 
ization of the (bosonic) Hubbard model [ij, which with 
nearest-neighbor hopping and onsite interaction is the 
minimal model of strong correlation physics 0] ■ 

Here, we describe and analyze numerically a particu- 
lar setup with ultracold bosons of two species in an opti- 
cal superlattice, described by a Bose-Hubbard model. In 
the limit of strong onsite interactions, the system can be 
described by the spin-1/2 Heisenberg antiferro- or ferro- 
magnet, depending on the parameters of the superlattice. 
The motivation is fourfold: 

(i) In the vein of Feynman's idea to simulate quantum 
systems by other quantum systems Q , it would be a great 
achievement to implement models of magnets like the 
Heisenberg model with ultracold atoms in optical lattices. 
In condensed matter systems, collective magnetism arises 
from the Coulomb interaction and the particle statistics 
which cause (super)exchange processes [1, H, @, 0, B H, 



[Toj . In particular, exchange interactions resulting from 
second order hopping processes in the Fermi-Hubbard 
model dominate its behavior in the Hmit of strong onsite 
interaction and are captured by an effective spin model, 
namely the Heisenberg antiferromagnet El, ■ 

While collective magnetism has been widely studied 
in solids over the decades, several experimental restric- 
tions apply quite generally: It is generally far from clear 
to what extent the typical simplified models are quanti- 
tatively realistic, and how to obtain the interaction pa- 
rameters. External control of these parameters is very 
difficult. Moreover, quantum magnetism becomes par- 
ticularly interesting in low dimensions. In real effectively 
low-dimensional solids it is however hard to control or 
assert the effect of the weaker interactions in the second 
and/or third dimension. Last but not least, solids give us 
only access to the linear response regime as sampled e.g. 
by neutron scattering. Questions of out-of-equilibrium 
many-body dynamics are essentially inaccessible. 

Experiments with ultracold atoms in optical lattices 
constitute clean and well-tunable manifestations of the 
Bose- or Fermi-Hubbard model [l|, [l^. To implement 
magnetic systems, the most straight forward approach 
would hence be to use a gas of ultracold fermions. How- 
ever, cooling of fermionic gases to the quantum regime 
is a considerably harder task due to the lack of s-wave 
scattering among identical fermions [ij, [l^ [l^, [13, [3 ■ 
Hence it is desirable to develop an alternative route via 
gases of ultracold bosons. Our investigations follow this 
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idea [H, HO, [HI, [H, [111 . Although we focus here on one 
dimension, analogous setups in higher dimensions 
could be used to investigate a plethora of frustrated spin 
systems that are hard to access analytically and numer- 
ically. 

(ii) The superlattice structure chosen in our setup (in 
analogy to the recent experiments [13, [IB]) allows in con- 
trast to [13, [2l[ for the tuning of the effective spin-spin 
interaction by changing an alternating scalar potential A; 
Fig [1] With the hopping strength t and the onsitc inter- 
action U of the Hubbard model, the coupling in the cor- 
responding effective spin model is then 4t^[//(L/^ — A^). 
This allows on the one hand to switch for the effective 
model between the Heisenberg ferro- and antifcrromag- 
net. On the other hand one might hope to increase for a 
fixed onsite interaction U the effective coupling, by choos- 
ing A ^ U. In this case, the relevant physics would 
become visible at correspondingly higher temperatures. 
However, the validity of the effective model breaks down 
in the vicinity oi A U. So one has to balance the va- 
lidity of the Heisenberg description and the temperatures 
needed to observe the quantum effects. To this purpose, 
the parameter A can be easily varied and used to tune 
to the Heisenberg regime in controlled fashion. 

(iii) In recent experiments [l^, [1^ by the Bloch group, 
the same optical superlattice as the one discussed here 
was used. But its parameters were chosen such that 
the superlattice decomposed actually into isolated dou- 
ble wells. The experiments analyzed dynamics in these 
double wells and contrasted in particular first order (hop- 
ping) processes (Hubbard regime) versus slower second 
order processes (Heisenberg regime). The next step 
would be to observe coherent many-particle dynamics af- 
ter coupling the double wells. We analyze such a situa- 
tion by the time-dependent density matrix renormaliza- 
tion group method (DMRG) [13, [2il. We focus on the 
coherent evolution of an initial Neel state and the differ- 
ences between the Heisenberg and the Hubbard regimes 
and present the experimentally available observables. 

(iv) Besides testing the coherence in the experiments, 
the setup allows to address questions of non-equilibrium 
many-particle systems, which is in general difficult for 
all present analytical and numerical methods. Contrary 
to the setup of isolated double wells, one observes for 
the many-particle dynamics in our setup a relaxation of 
local quantities. This is an indicator for convergence 
of subsystems with finite real-space extent to a steady 
state. Recently, the mechanism of how such a relax- 
ation may occur was clarified for (free) integrablc sys- 
tems in |29| . Corresponding examples can also be found 
in [13, l3ll. |32|. [33} . For a few nonintegrable systems the 
question was anal yze d numerically in |33 . [ssl . [36j and an- 
alytically e.g. in |37l . [s^ . In general one expects that 
in nonintegrable models, thermalization occurs (due to 
scattering processes), and that in integrable models, re- 
laxation (to a nonthermal steady state) occurs via phase 
averaging effects [1^. This is demonstrated here ana- 
lytically for the Heisenberg model. Our setup could be 



used to study experimentally such relaxation processes 
- in particular, the qualitative differences between non- 
integrable systems, here the Bose-Hubbard model, and 
Bethe ansatz integrable models [13, [13] , here the Heisen- 
berg model. 

We also fill a certain gap of current literature on 
such topics (see e.g. [13, [Hi, [HI) by emphasizing that 
the Heisenberg spins of the effective model, obtained 
by the Schrieffer- Wolff transformation (in Appendix |^ 
[41| . should not be identified directly with the two bo- 
son species. A spin up of the effective model corresponds 
rather to a particle of species I dressed by hole-double- 
occupancy fluctuations. The analogy holds only for small 
t/{U ± A). The consequences for experimentally avail- 
able observables are surprisingly strong. Recently in [13| 
a setup of coupled double wells (A = 0, but alternating 
hopping t 7^ t' in Fig. [1]) was analyzed numerically - in 
particular, the possibilities to generate entangled pairs 
of particles were studied. Again, a perfect mapping to a 
spin model was assumed from the outset. The validity of 
this mapping is one central topic of our article. 

The paper is organized as follows. Section [II] de- 
scribes the experimental setup and how it can be de- 
scribed by a Bose-Hubbard model. Restricting to half 
filling. Section Hm and Appendix [K\ derive by a Schricffer- 
Wolff transformation for the limit of large onsitc in- 
teractions an effective model which is the Heisenberg 
antiferro- or ferromagnet. In Section llVl we investigate 
by time-dependent density-matrix rcnormalization-group 
(DMRG) the evolution of typical observables like mag- 
netization, momentum-space and real-space correlators, 
where the first two are also available experimentally. The 
focus is on contrasting the differences between the full 
Hubbard dynamics and the corresponding effective spin 
model, and also the differences to the case of isolated 
double wells [13, [13] • We observe indications for (lo- 
cal) relaxation to steady states. This is discussed in Sec- 
tion |V] where we also explain how the relaxation for the 
Heisenberg model is connected to a phase averaging ef- 
fect. Section IVII addresses in more detail the question 
why and under what circumstances the effective model 
is a valid description for the full Hubbard Hamiltonian, 
especially for the dynamics. In Section |VlT] we argue that 
the groundstate of the Heisenberg antifcrromagnet could 
be prepared by tuning an alternating hopping parame- 
ter of the superlattice adiabatically. Section IVlIII gives a 
short conclusion. 



II. SETUP AND MODEL 

In the following, we present a setup of ultracold bosonic 
atoms in a one-dimensional optical superlattice that re- 
duces in certain parameter regimes, where first order hop- 
ping processes are suppressed, to the Heisenberg ferro- 
or antifcrromagnet. The use of bosons is motivated by 
the fact that experimentally, access to the low energy 
quantum physics is at the moment still much harder for 
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FIG. 1: a) Tight-binding system parameters for our system 
of ultracold bosons of two species (t, I) in a one- dimensional 
optical superlattice ([T}. b) If the phase shift between the 
two laser potentials (dashed lines) is zero, the minima of the 
full potential (solid line) are equidistant, the Wannier wave 
functions for each site are reflection symmetric, and hopping 
parameters t and t' are hence equal. 



fermionic systems. In analogy to the fcrmionic case, for 
which the antiferromagnetic Heiscnberg model describes 
the effective low-energy physics of the fermionic Hubbard 
model, we choose to have two species a G {t, i} of bosons 
in the lattice - two hyperfine states of a bosonic atom. At 
half filling {N sites, = Ni = N/2), the effective low 
energy model is, as we will see later, the ferromagnetic 
Hcisenberg model. To allow for tuning the effect of first 
order processes and to switch between a ferromagnetic 
and antiferromagnetic regime, we employ an alternating 
onsite potential A,; and call the two sublattices A and B. 
The potential minima differ by a value A > 0. Such su- 
perlattices can be generated by the superposition of two 
laser frequencies of ratio 1:2, see Ref. [lH, H^, \^ and 
Fig. [TJ Further, the tight-binding approximation with 
restriction to the first Bloch band (one Wannier function 
per site) is assumed. Then the system is described by 
the two-species Bose-Hubbard Hamiltonian 



^ a 



/i.e.) +J2^' 



U^nT,ni, + ^^naiinai~ 1), (1) 



In particular, we choose 

^ I —A/2 for even i 
' ^ I A/2 for odd i. 



The superlattice potential is of the form 

V{x) = Vo ain^ikx) + Vi sin^{kx/2 + , 



(3) 



The amplitude Vi of the second potential can be used to 
tune A; for our purposes, Vi ^ Vq. In ([1]) it was as- 
sumed that hopping from a site i to its neighbors i ± 1 
occurs with equal amplitude t. In principle, the hop- 
ping depends exponentially on the distance between the 
potential minima [l], [4^. Only if we choose the phase 
difference to vanish as in Fig. [T|3, the positions of the 
potential minima will be equidistant for all Vi, and the 
Wannier wave functions for all sites are reflection sym- 
metric. In this case, the hopping will hence be of equal 
strength for all bonds, i.e. t = t' in Fig. [T] This situation 
is considered in the following. 

For the actual analysis of dynamics in Section IIVI 
we will choose equal inter- and intra-species interaction, 
U ~ Ug. This is at the moment the standard situation 
for the corresponding experiments. With U = Us, the 
effective Heisenberg models, describing the second order 
physics, will turn out to be isotropic. The onsite poten- 
tial A and interaction U can be calculated in harmonic 
approximation of the lattice potential around its minima. 
For = in ([3]) , the corresponding oscillator frequencies 
are nw± = 2^Er{yo ±Vi/A) « 2/^^(1 ±i^), where 

1-2 7 2 

Er = -^^^ is the recoil energy of the laser potential with 
the shorter wave length. For the dependence of the onsite 
potential on the lattice parameters follows 



A = V^i-f(^+-^„)«yi 



(4) 



In order to achieve an effectively one-dimcnsional lattice, 
V(x) is superimposed with two transversal laser beam po- 
tentials of a (higher) amplitude V±_ . Within the harmonic 
approximation, and with the s-wave scattering length a^, 
the resulting onsite interactions for two neighboring sites 
are El 



(K) ± Vi/A)Vl 



(5) 



(2) 



They are in principle not identical. Irrespective of this, 
the effective spin model derived in Section Hill would be 
isotropic and translation invariant. As we will see in the 
following for a set of realistic experimental parameters, 
we have for the case Vi = (=> A = 0) that Vq ^ Er 
and Vb ^ U . The onsite potentials A considered in our 
numerical simulations of the Hubbard model ([1]) are from 
the interval A G [0,4[/]. According to (g]), Vi will for 
nonzero A hence obey Vq S> and we can thus use 
U+ = U- = [/ in good approximation. 

For the tight-binding approximation to hold, we need 
that energies t, U , and A are well below hui± « 2-\/E7^o, 
the energy scale for vibrations of an atom in one min- 
imum of the laser potential. With the hopping t = 

^ {VoVl/E^rf'\-^V^^Er [il,[43, we have for ex- 
ample with A = 27r/fc = SOOnm, Rubidium atoms (i.e. 
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« ttO.OI), ^0 « 8.7Er and V± « 30^;^ (cmp. e.g. to 
[111) that t ^ OMEr-, 17 « 8t and Huq ^ lOOi. So as long 
as A is also well below lOOt, the tight-binding approxi- 
mation for ^ using the lowest Bloch band is valid. This 
will be the case in the rest of the paper. 

III. EFFECTIVE MODEL 

To go to a regime where the physics of the two-species 
Bose-Hubbard model ^ reduces to that of a Heisenberg 
magnet, we choose half filling 

N^^N^= N/2, (6) 

{N is the total number of lattice sites) and assume the 
large- C/ limit 

t<|C/±A|. (7) 

In this limit, occupation of a single site by more than one 
boson is energetically unfavorable and will occur only in 
short lived intermediate states. This means that single 
(first order) hopping processes are suppressed. Besides 
some hybridization effects, exactly one boson sits on each 
lattice site and can be identified with an effective spin on 
that site (up and down orientations corresponding each 
to one of the boson species). The second order hopping 
processes as depicted in Fig.[5]lead then to nearest neigh- 
bor spin-spin interactions. 

We will show below, how the effective Hamiltonian can 
be derived by a Schricffer- Wolff transformation. While 
this is a well-known procedure, there is an interesting 
twist to the interpretation of the result. For the mo- 
ment let us work with the naive identification of spins 
up and down of the effective model with the two boson 
species of the full model ^ . We want to derive an effec- 
tive Hamiltonian describing the physics of the Hubbard 
Hamiltonian ^ in the subspace Tii of singly-occupied 
sites. 

7^1 :=span{|T>,U>r^ (8) 

The effective Hamiltonian can be deduced from the fol- 
lowing simple recipe: With exactly one spin per site, the 
onsite interaction is ineffective. Hopping processes occur 
only in second order, leading to a spin-spin interaction 

The corresponding coupling strengths J and Js are ob- 
tained by dividing for each possible second order process 
(Fig. ^ the product of the transition matrix elements t • t 
for the hopping to a neighboring site and back by the en- 
ergy difference [/ ± A (Us ± A) to the intermediate state, 
and adding all such terms that contribute to the same 
effective spin-spin interaction (see also H^). 

WTa ^ Tj~~K ^ c/2 - a2 ' '^-^uf~^- 

(10) 



The effective Hamiltonian ([9]) is the XXZ model. From 
now on we specialize to U = Us, i.e. J—Js = — J and have 
hence the isotropic Heisenberg ferromagnet for A < U 
(J > 0) and the isotropic antiferromagnet for A > U 
(J < 0). Note that the effective Heisenberg Hamiltonian 
would also be isotropic and translation invariant if the 
onsite interaction U would be different for even and odd 
sites. If the hopping would be alternating (t 7^ t'), we 
would obtain the dimerized Heisenberg model. 
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FIG. 2: Some second order hopping processes (superex- 
change) contributing to the effective spin model Q when first 
order hopping processes in the full Hubbard model ([1} are 
suppressed. The figure displays possible initial states (light 
color) with one particle per site and intermediate states (dark 
color) with doubly occupied and empty sites. 

A common mathematical approach for the deduction 
of such effective models is the Schricffer- Wolff transfor- 
mation [111 . We are interested in the physics of the sub- 
space Til with exactly one particle per site. The Hub- 
bard Hamiltonian couples this subspace in first order of 
the hopping t to the rest of the Hilbert space (states with 
doubly occupied and empty sites). The Schricffer- Wolff 
transformation, 

HX ■■= e'^He'^^, nf^ e~'^ni, (11) 

is a unitary transformation with generator S chosen such 
that the transformed Hamiltonian does not contain 
terms anymore that couple Tii to the rest of the Hilbert 
space, or at least only in some higher order of t. In 
Appendix lAl a generator 

is derived, such that effective Hamiltonian is 

iv) 

The full effective Hamiltonian (|A16p still contains a 
term i[S, H^] representing the remaining coupling of the 
subspace Hi to the rest of the Hilbert space which is of 
order t^. 
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The method is based on the smallncss of S. Accord- 
ing to (|12p . it hence breaks down when U ~ |A|. The 
effective Hamiltonian ([T3|) is only vahd for |C/ ± A| ^ I. 
Only in this regime, the first order hopping processes 
leading out of Tii are suppressed. See the discussion in 
Section |VD 

Often spin up (down) states of the effective model are 
then identified with a boson of species t (species |) on 
the corresponding sites. However, with respect to the 
original model, Eq. ([1]), it is not Hi itself that is weakly 
coupled to the rest of the Hilbert space and evolves ac- 
cording to the Heisenbcrg Hamiltonian but the subspace 
H™'^ defined in Eq. pT|). A spin up in the effective model 
corresponds rather to a boson of species t with a cloud 
of hole-double-occupancy fluctuations a^i — > e^^ aaie~'^^ ] 
see also Fig. |31 The experimental consequences arc sur- 
prisingly strong, as we will in the next section. 

IV. TIME-EVOLUTION FROM THE NEEL 
STATE 

In the following, numerical results for the evolution of 
the system where the initial state is the Neel state 

I'/') := I Tititi • . ■> e (14) 

are presented. This parallels recent experimental inves- 
tigations [1^, [2^ of the evolution of corresponding states 
I ti) in isolated double wells. 

To compare the effect of first and second order pro- 
cesses, the evolution was done twice for each set of pa- 
rameters {i = 1, U = 8, various A; see Section [Til) , oi^ce 
with the corresponding Heisenberg model (in the sub- 
space Til) and once with the full Hubbard Hamiltonian 
(in the full Hilbert space), where the initial state (fH)) 
was in fact chosen as the tensor product of alternatingly 
having one boson of species t or j on each site. The two 
different time scales of first and second order processes 
become clearly visible. The qualitative differences to the 
isolated double well situation (as analyzed in [2^, [2^) 
and resulting interesting questions for experimental in- 
vestigations are discussed. 



A. Errors through experimental limitations in 
state preparation and measurement 

We shortly want to discuss how well the dynamics 
of the magnetic model, the Heisenbcrg model, can be 
implemented experimentally by those of the two-species 
Bose-Hubbard model. In the literature on magnetism 
via ultracold two-species atom gases in optical lattices 
[20l . [2ll . [23! , spins up and down of the magnetic system 
are usually identified directly with atoms of species t 
and [ of the ultracold gas. In this vein, evolution of the 
Neel state with the Heisenbcrg Hamiltonian would 
be translated into evolution of the state | tititi • • •) 



with the Hubbard Hamiltonian. This is actually correct 
only to zeroth order in S . 

We want to implement the evolution of a state 
\(f)) £ Hi, under the effective Hamiltonian Hes = 
e^'^He~'''^\-Hi by the evolution of a state Itp) = e~^^\(l)) € 
T^i"^, (fTT|) . under the Bose-Hubbard Hamiltonian H. 
The state \^) is \(j)), superimposed with states where 
starting from pairs of doubly occupied sites and 
empty sites were created (e.g. Eq. and Fig. [3] be- 
low). This can also be interpreted as constructing the 
Neel state with effective spins, each corresponding to a 
boson accompanied by a cloud of hole-double-occupancy 
fluctuations a^i — *■ e^^ a^iC^^^ . The decisive point is now 
that it seems not possible to prepare such states from 
7i"'® (and has to our knowledge never been done), but 
only some specific states from Hi- Hence, instead of 
starting the experiment from the initial state \'4)), one is 
forced to neglect the Schrieffer- Wolff transformation and 
start from the state - in our example the Neel state. 
For observables that do not change the number of doubly 
occupied sites, this results in an error of 0(5^). 

If we had determined the exact Schrieffer- Wolff trans- 
formation e"'"^ (i.e. S exact to all orders in t) and could 
actually implement it e.g. by time evolution in the ex- 
periment, all measurements would be exact. One could 
prepare the state |</>), apply the Schrieffer- Wolff transfor- 
mation by time evolution to obtain evolve with the 
Hubbard Hamiltonian for some time i, apply the inverse 
Schrieffer- Wolff transformation and measure our observ- 
able O. This would yield the exact equality 

((^|e-^<=«*/'''-0-e^<^«*/'''|(/)) 

= {4>\e'^e-"'''^e-'^ ■ 6 ■ e'^ e"'''^er'^\4>), (15) 

where H^s would now of course be a generalization of the 
Heisenberg model with longer ranged interactions. 

In the Appendix I A 11 S is determined to first order in 
jj^^ , (|A15[) . and correspondingly i?off to first order in 
the effective coupling J, (|A21|) . Using this approxima- 
tion of iS instead of the exact one, the remaining errors 
in the observables are of order S'^ , i.e. (It is 

not 0(5^), because the operator 5, given in Eq. ([Alii) , 
changes the number of double occupancies by one and the 
typical observables O we are interested in do not.) How- 
ever, failing to implement the Schrieffer- Wolff transfor- 
mation completely, i.e. measuring ((/)|e~^*/'''-0-e^*/'''|(^) 
instead of p5|) . leads to errors of order . This will be 
demonstrated in an example (Section IIVDI Figs. [3 and 
[HI). In addition to the error from neglecting or truncat- 
ing the Schrieffer- Wolff transformation, there is the error 
from truncating the effective Hamiltonian (|A2ip . This 
accumulates with time and is in principle of order ,Pt, 
but may also just result in a sort of rescaling of the time 
axis. Also the local observables considered relax rela- 
tively quickly, making this second source of error less 
important. 
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In the remainder of the article, the initial state ([T4|) . 
evolved with the Heisenberg Hamiltonian will be 
called (j){t). If it is evolved with the Hubbard Hamil- 
tonian it will be called (p{t). 

:= e^-*/»'"'|0) and |0(O> := e^*/^''|0), (16) 
and concerning observables we have explained that 

To illustrate the considerations above, let us shortly 
regard the case of an isolated double well (two sites). 
Hamiltonian and S read (cf . Appendix |^ 

H = Ht + Ho 

Ho = (C/ + A)|0,ti)(0, ti \ + iU-A)\ ti,0)(ti,0| 

^t = -t(|o,Ti>(T,il + ITi,o)(i,t| + /^.c.) 

'5-i7^in,0)((i,T| + (T4l) 

-u^\OAi){{UU + {tM)+h.c.. 

With this, the effective Hamiltonian and the transformed 
initial state are 

^^cff = —j^z^^i ■ ^2 + 0(5"*) 

=^iTa) + T7^in,o)-^io,n). (is) 

So a magnetic state with one particle per site corresponds 
in the experimentally realized Hubbard model to the 
magnetic state plus an admixture of states with doubly 
occupied and empty sites, Fig. [31 The original Hamil- 
tonian generates with Hi doubly occupied sites to first 
order in t. Conversely, in the (full) effective model, such 
terms are at least of order (in the two site case here, 
actually of order f*). 



Heisenberg | | two species Bose-Hubbard 




FIG. 3; Spin up (down) states of the effective magnetic 
model are not to be identified directly with a boson of species 
t (species |) in the experimentally realized Bose-Hubbard 
model. A spin up in the eff'ective model corresponds rather 
to a boson of species t with a cloud of hole-double-occupancy 
fluctuations aa-i — > e^'^flo-ie^''^. In the vicinity of A = [/, the 
correspondence breaks down. 



B. Symmetry between the ferromagnetic and the 
antiferromagnetic cases 

The Neel state (dH), and the effective Hamiltonian 
(fO)) are both real in the {5f }i-cigenbasis B :— {\cr) = 
\<Tia2 ■ ■ •)} (real coefficients and matrix elements). Typ- 
ical observables O of interest like for the magneti- 
zation or SfSi and S^S~ + StS~ for correlators are 
real in that basis and selfadjoint. It follows that the cor- 
responding expectation values {0)^i^t) are identical for 
the Heisenberg ferromagnet (J = 1) and antifcrromag- 
net (J = —1): Let Oo-^ct' := {(t\0\(t'), := {(T\(f)), and 

M(a)fm,<T.<T'(0 := (cr|^^(a)fm(i)k') for the time evolution 
operator of the (anti)ferromagnetic Heisenberg model. 
Then 

= 0^(4nW)*-O-KmW)*'/' 

= mltjt)-o-u,ijtM)- (19) 

The evolution of the corresponding observable under 
the full Hubbard Hamiltonian {O) ^^^^ will obey this sym- 
metry to zeroth order in S. Typically, the resulting 
curve will coincide well with the corresponding Heisen- 
berg curve. The smaller jC/^ — A^j is chosen, the worse 
the effective model will capture the actual dynamics and 
the stronger deviations from the corresponding Heisen- 
berg results will be. The specific form of the devia- 
tions, however, will depend on the choice of U , A, and 
t. In particular they show no symmetry when switch- 
ing between the antiferromagnetic and the ferromagnetic 
regimes (A ^ U). To illustrate this further, several plots 
contain the two curves A = and A = y/2U which have 
according to (jlO[) the same effective spin spin interac- 
tion strength J, except for the opposite sign (FM, AFM, 
respectively). 

C. Numerical method and parameters 

For the numerical simulation, a Krylov subspace vari- 
ant |44|. [isl o f the time-dependent DMRG algorithm was 
used [27l. l28l. l46j. For the Hubbard model, the site basis 
was restricted to a maximum of two particles for each 
species. Insensitivity of observables to the chosen maxi- 
mum number of bosons per site was affirmed. We chose 
lattice sizes of i = 33 for the Bose-Hubbard model and 
L = 65 sites for the Heisenberg model. Odd numbers are 
useful here to have reflection symmetric states. 

As a matter of fact, boundary effects are much less 
problematic here than in groundstate calculations, as the 
initial state is a product state and correlations between 
sites are generated inside a causal cone (the analogon 
of a light cone; see Section IIVEI in particular Figs. [TU] 
and I12p . So as long as measurements are done in the 
middle of the system, outside of the causal cones starting 
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from the boundary sites, results are except for exponen- 
tially small contributions identical to those of an infinite 
system (thermodynamic limit). 

In the time evolution, the absolute difference per phys- 
ical time unit between exactly evolved state and the state 
evolved via DMRG 1 1 V'Sr''* " i'dt^^^ I |/diiV was bounded 
from above by e = 10~^ to £ = 10~^ and the time step 
chosen appropriately between dt — 0.1 and 0.01. The 
errors were determined in a rigorous fashion, by calcu- 
lating the exact value of || |fc -I- 1) — H\k) ||, where \k) 
are the Krylov vectors. For all calculated observables, 
convergence in the error bound and dt was checked. The 
resulting number of basis states, used to represent the 
time-evolved state, was < 3000. 



D. Site magnetization 

Figures I3H51 show the evolution of the site magnetiza- 
tion nix ~ {Sx)(t> in the Heisenberg model and the corre- 
sponding quantity m^=(n|^-n|^)^/2 (= (S^)^+0{S'^) 
according to Section |IVA[) for the full Hubbard Hamilto- 
nian. For the latter, times were rescaled by the cou- 
pling constant J, p^ . of the corresponding effective 
spin model. Site x was chosen to be in the middle of 
the system in order to avoid finite size effects (cf. Sec- 
tion [IVC]). For the Heisenberg model (in the thermody- 
namic limit), the site magnetization obeys for symmetry 
reasons mx+i = —irix- Analogously, due to invariance 
under translations by an even number of sites and parti- 
cle number conservation, one has for the Hubbard model 
(again in the thermodynamic limit) (n^x+n-a-x+i) ^ = IVt 
and hence rhx+i = -rhx, and {n^x+i + '^ix+i)^ = 
2 — {n^x + "i^lx)^ for all times. As discussed in Sec- 
tion IIV C[ deviations of our numerical results from the 
thermodynamic limit arc negligible, although the simu- 
lations are carried out with finite lattice sizes. 

The larger jC/^ - A^] is (for fixed t = 1), the bet- 
ter the curves for the full Hubbard Hamiltonian coin- 
cide with those of the Heisenberg model. This is consis- 
tent with Sec. HIT] as the perturbative derivation of the 
effective model becomes exact in this limit. Note that 
the deviations between the measurements stem from two 
contributions here: (a) failure of preparing the correct 
7i"'® state, i.e. applying the Schrieffer- Wolff transfor- 
mation at t = 0, (b) failure of measuring 5*1 instead of 
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S"! - i[S, 5|] + 0(52). The weight of those 
errors which are of order vanishes only far away from 
|A| = [/. 

For the Heisenberg model, we observe relaxation of the 
site magnetization from ±1/2 to 0. The oscillations of 
this observable occur on the time scale 1/J. The relax- 
ation is possible due to the continuous spectrum of the 
Heisenberg model (in the thermodynamic limit). Here, 
the convergence to a steady state is connected to a phase 
averaging effect, as is typical for integrable systems. An- 
alytically this can be seen in a time-dependent mean 



field treatment of the Heisenberg model which we present 
in Section IV Bl For the staggered magnetization [mx) 
one obtains (in this approximation) a damped oscillation 
with the amplitude decaying as ~ l/t^/"^. This coincides 
well with the DMRG data, giving support to the mean 
field approach; see Section IVBl Fig. [T51 

For large \U'^ — A^| the Hubbard dynamics clearly fol- 
low the curves obtained with the Heisenberg model (sec- 
ond order processes); Fig. |4]and[5l On the shorter time- 
scale 1/t = 1 (J/t in the rescaled plots), correspond- 
ing to first order processes, small oscillations around the 
Heisenberg curves are visible. Their amplitude decreases 
with increasing |[/^ — A^|. The perturbative treatment 
of the system, leading to the isotropic Heisenberg model, 
breaks down for |A| ~ [/. In this case, the two boson 
species cannot be interpreted as spin up or down states 
anymore and one has an appreciable amount of double 
occupancies in the system as demonstrated in Fig. [6l 



Hubbard, A = 8 V2 
Hubbard, A = 16 
Hubbard, A = 24 
Hubbard, A = 32 
AFM Heisenberg 




FIG. 4: Evolution of the magnetization on a particular site x, 
starting from the Neel state and evolving with respect to the 
full Hubbard Hamiltonian with i7 = 8 and A > [/, and the 
isotropic Heisenberg antiferromagnet, respectively. The first 
order processes occur on the time scale t = 1 (here 1/J due to 
the rescaling of the time axis, where time is given in units of 
the effective coupling J) and their amplitude decreases quickly 
with increasing |f/^ — A^|. 

Finally, we want to compare those results to the dy- 
namics for isolated double wells as addressed experimen- 
tally in [2^, [2^ . Fig. [7] shows for this case the dynamics 
of the site magnetization again for the Hubbard model 
at various A > C/ and the corresponding antiferromag- 
netic Heisenberg model. The decisive difference is that 
no equilibration is possible in this case. This is due to the 
fact that the Hamiltonian has only a few discrete eigen- 
values here, as opposed to a gapless continuous spectrum 
for the lattice systems in the thermodynamic limit. In 
the two-site Heisenberg model we have only two states 
in the basis of the 5^ = Hilbert space. The two eigen- 
states have energy difference J. The magnetization curve 
for the Heisenberg curve is hence just a cosine with fre- 
quency J and constant amplitude 1 . The dynamics of the 
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Hubbard, A = 
Hubbard, A = 2 
Hubbard, A = 6 
FM Heisenberg 
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FIG. 5: Evolution of the magnetization on a particular site 
X, starting from the Neel state and evolving with respect to 
the full Hubbard Hamiltonian with U = 8 and A < {/, and 
the isotropic Heisenberg ferromagnet, respectively. The first 
order processes occur on the time scale t = 1. Here the contri- 
butions of the first order processes cannot be made arbitrarily 
small as we are limited by |f/^ - A^l < The Heisenberg 
curve here is identical to the one of the antiferromagnet in 
Fig. |4] due to symmetry, see Section IIVBI The effective cou- 
pling J has the same modulus for A = and A = \/2U , 
namely | J| = 4t^/C/, but opposite sign. The two curves show 
quite different behavior. There is no particular symmetry 
except the one for the second order physics as discussed in 
Section HVB] 



corresponding two-site Hubbard model is determined by 
three (discrete) incommensurate frequencies. The mag- 
netization is hence not completely periodic, but due to 
the relation to Heisenberg model, a frequency ^ J is still 
dominating. No sign of relaxation is visible. 

As mentioned above, the differences between Hubbard 
and Heisenberg dynamics stem from the fact that the 
two Schrieffer- Wolff transformations in pV Ap have been 
neglected. Fig. [7] shows in the lower panel the site mag- 
netization for the case where both error sources (a) and 
(b) have been corrected. Although this should be hard 
to implement experimentally, it is unproblematic for our 
numerical analysis. We apply the Schrieffer- Wolff trans- 
formation (jA14[) . correct up to C((7^^), to the initial 
state before the Hubbard time evolution and its inverse 
before the measurement. As discussed in Section IIV A| 
the remaining deviations from the Heisenberg curve are 
then only of order 5"^; Figs. [7] and [51 



Correlation functions 



The correlation functions in Figs. support on the 
one hand the results already obtained from the magne- 
tization dynamics in Section IIVDI On the other hand 
one also sees here explicitly that correlations spread out 
inside a causal cone (analogon of a light cone) defined 
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FIG. 6: Evolution of the occupation number {n^x + nix) — 1 
(upper panel) and its variance (n^) — (ux)^ (lower panel) on 
a particular site x, starting from the Neel state and evolving 
with respect to the full Hubbard Hamiltonian with U = S 
and several A. The two quantities should be exactly zero, if 
the analogy to the spin model was exact. The analogy breaks 
when |[/^ — A^l goes to zero. In the special case A = 0, the 
system is (additionally to the invariance under translations 
by two sites) invariant under translation by one site plus in- 
terchange of particle species. Hence {n-jx + riix) = f Vt for 
A = 0. 



by the maximum group velocity. The latter coincides for 
large |C/2 - A^j with the maximum group velocity 2 J of 
the Heisenberg model. One also notes here that equilibra- 
tion to a steady state occurs first for small subsystems. 
This issue will be discussed in Section IVl 



F. Momentum distribution and correlators 

Experimental access to onsite magnetization or, corre- 
spondingly, the particle number difference (^pin imbal- 
ance" ) has already been demonstrated [1^, [2^ . However 
there is no direct access to the real-space correlators. As 
it turns out, the standard experimental observable for 
experiments with ultracold atoms, the momentum distri- 
bution (rife) = {nijk+nik) , is to zeroth order in S constant 
in time, ft measures to this order simply the particle den- 
sity which is in the limit of Heisenberg dynamics A 3> t/ 
very close to one. 
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(20) 



It follows with nil) 

{nk)^^{nk)4, + 0{S')^l + 0{S'). (21) 
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FIG. 7: Evolution of the magnetization (ttii := (n|i — 
rij i)/2) on one site of an isolated double well, starting 
from the Neel state and evolving with respect to the full 
two-site Hubbard Hamiltonian with U = 8 and A = 
10,16,24,32, and the isotropic Heisenberg antiferromagnet, 
respectively. Contrary to the case of an infinite lattice, 
the magnetization does not relax here. The upper panel 
shows (mi)^ = {(l)\e-"^^'''mie"^^'''\cl>). In the lower panel 
shows {(j)\e'^e-"'^'''e~'^mie'^e'^'^'''e-'^\(l)}, i.e. there the 
Schrieffer- Wolff transformation was accounted for (5 correct 
to 0{jj^-^)). As discussed in Section llV Al the stretching 
in the curves w.r.t. time results from terms of order in 
the effective Hamiltonian. They originate from fourth-order 
hopping processes. 
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FIG. 8: Time average of the particle number variance on one 
site of an isolated double well (ni = n|i +nii), evolving with 
respect to the full two-site Hubbard Hamiltonian with several 
U and A (J = 4t^U/{U'^ - A^)). The variance should be 
exactly zero, if the analogy to the spin model was exact as we 
would have exactly one particle per site then. As discussed 
in Section HVH the error is of 0(5^), if the Schrieffer-Wolff 
transformation is neglected completely, (|17|) . and of 0(5*), if 
its first order approximation (|A14p is used. In the special case 
of the isolated double well, the second order terms in 5 vanish 
(because Ht = here, see (|A6|l ). Hence, we actually observe 
0{S^) instead of 0(5"). The quantity on the a;-axis quantifies 
0(5). For each curve, either A ot U was kept constant and 
the other parameter varied. Compare also to Figs. 3] and [S] 



Hence one needs to go beyond the measurement of 
the momentum distribution. By analysis of shot-noise 
of absorption images taken by time-of-flight measure- 
ments, one obtains momentum-space particle density cor- 
relation functions [13, . Experimentally available arc 
{nakUak') and {nkUk') = Y.aa'i'^'^k'ricr'k') and hence also 
X^o-i'^o-fc'^-CTfc')- In the following, we will again use the 
approximation (6)^ = {6)^ + 0{S^) ~ (O)^. 

(aLa<Tjal'mVn)0 (22) 

Those give information about long-range spin correla- 
tions. 

+ (1 - (5y)(5i„<5jma|ia|ja|ja^-)0 

= {5ij5mn{\ + S'f )(^ + S^^-i) 

+ (l-5,,)<5,„,5„„(i + ^f)(| + 5|))^ 

(23) 



for the observable {n^kn^y) and 

— {^ij^mn ^^(7 ^ai'^—am 

+ (1 - 5ij)5,n5jrn E<t "^aiO-ajO^-aja-ai) 4' 
= (%'5_(i - 2^f5,^) 

+ (1 - 5.,,)5„,5,r,^2{SfSf^ + SfSyj) ^ (24) 

for the observable '^^{n^k'n^ijk') ■ This also reflects the 
fact that to zeroth order of S, there are no double occu- 
pancies with respect to the original basis. However, there 
is an admixture of them, contributing in second order of 
S, (?iT»"i»>0 = 0(5^). Inserting ^ and ([24 ]) to ((22 l) 
yields 

{n^kn^k')^ -\ - -h + J^kk' + {Sl,k^^^Ak)<t> 

= {niknik')<p (25) 

and 

^{nakn-^k')^ 

(T 

— ^ ~ J! ^ '^(^Ak^-Ak + ^Ak^-Ak)4> (2^) 
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FIG. 9: Evolution of the analogon j{{nfx — n^x){'n'fx+e — 
nix+i))4, - K^Tx - nix)^{n^x+i - of the magneti- 

zation-magnetization correlation function, starting from the 
Neel state and evolving with respect to the full Hubbard 
Hamiltonian with (7 = 8 and A = 16. The plot shows the 
absolute value of the correlator in logarithmic scaling. The 
line denotes the maximum group velocity 2J of the Heisen- 
berg model. 
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FIG. 10: Evolution of the magnetization-magnetization cor- 
relation function {SxSx+e)tl>~{Sx)ri>{Sx+i)4>, starting from the 
Neel state and evolving with respect to the isotropic Heisen- 
berg antiferromagnet. The plot shows the absolute value of 
the correlator in logarithmic scaling. The line denotes the 
maximum group velocity 2J of the Heisenberg model. 



where Afe ^ k-k', 5^ ^ ^ e^-^-^^f , ^" ^ ^ 4", 
and S^\(f>{t)) ~ were used. 

A numerical comparison of the evolution of the 
momentum-space spin-spin (density-density) correlator 
for the Heisenberg and the Hubbard models is given in 
Figs. [TSlfTSl To achieve such a good agreement, two cor- 
rections were necessary that are described in more detail 
in Appendix [B] First of all one needs to correct for finite 
size effects. Secondly, single particle Green's functions 
(alaj) enter which are trivial when evolving with the 

Heisenberg model {{a\a^)^ = Sijn^i{t)), but have con- 
tributions of 0{S^), when evolving with the Hubbard 
Hamiltonian. In the comparison of both evolutions, they 
can hence be understood as a major carrier of distur- 
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FIG. 11: Evolution of short range magnetization- 
magnetization correlation function (distances i = 1,2,3), 
starting from the Neel state and evolving with respect to the 
full Hubbard Hamiltonian with (7 = 8 and A = 16 and the 
isotropic Heisenberg antiferromagnet, respectively. With in- 
creasing time, deviations between Heisenberg and Hubbard 
dynamics become more pronounced than for the magnetiza- 
tion in Fig. |4] However, in both cases tendency toward equi- 
libration to a steady state is visible. 
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FIG. 12: Evolution of the magnetization-magnetization cor- 
relation function, starting from the Neel state and evolving 
with respect to the full Hubbard Hamiltonian with (7 = 8 
and A = 16 and the isotropic Heisenberg antiferromagnet, 
respectively. See also Fig. 1111 



bance, reflecting first order processes in the Hubbard 
model. To achieve comparability it would be desirable to 
remove contributions from {a\a^) completely. This would 
be possible for our numerical analysis. In a correspond- 
ing experiment however, the quantities are not available. 
Hence we confined ourselves to removing only the contri- 
butions from nearest neighbor correlators {a\a^^■^. As 
Figures [T51fT5l demonstrate that this is already sufficient 
and Fig. [TC] that it is necessary. The experimental obser- 
vation of the nearest-neighbor correlators is within reach 
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The specific form of the momentum-space correlation 
function can be understood with the causal cone behav- 
ior of the corresponding real-space correlators discussed 
in Section [IV El At the beginning of time evolution, cor- 
relations for small distances build up (e.g. due to the 
spin flip terms S^^.-^ in the Heisenberg model). This 
corresponds in the momentum space representation to 
correlations for large Ak. As the correlations spread out 
in real-space, correlations for smaller momenta Ak build 
up. 
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FIG. 13: Evolution of the momentum-space spin-spin corre- 
lator {5'Afc5'lAfe)0 for the Heisenberg antifcrromagnet. The 
correlator corresponds according to Eq. p5|) to the density- 
density correlator in the Hubbard model and is available in 
experiments with ultracold atoms [47| . The initial state (|14p 
is uncorrelated. Correlations build up on the time scale 1/J. 
Finite-size effects have been corrected (see text). 



V. RELAXATION TO STEADY STATES 

A. General features 

Contrary to the setup of isolated double wells, one ob- 
serves for the many-particle dynamics in our setup a re- 
laxation for local quantities. This may be seen as an in- 
dicator for convergence of the states of subsystems with 
finite real-space extent to a steady state. Recently, the 
mechanism of how such a relaxation may occur was clar- 
ified for (free) integrable svstems | 2 9| . Corresponding 
examples can also be found in |30l. l3ll. Is^ . [33[. The 
setup considered in this paper could be used to study 
experimentally such relaxation processes - in particu- 
lar, the differences for the nonintegrable Bose-Hubbard 
model and the Bethe ansatz integrable Heisenberg model. 
Experimental investigations would be very useful here, 
as the fast entanglement growth during time evolution. 
Fig. 1171 prohibits numerical access to long times and for 
Bethe ansatz integrable systems, analytical results arc 
also relatively limited for such purposes. 




FIG. 14: Evolution of the momentum space density-density 
correlator '^f^irifk+Akni^k) ^ (minus the trivial parts on 
the right hand side of Eq. (|25|) ) in the Hubbard model with 
?7 = 8 and A = 16. Except for quick oscillations on the 
time scale t — 1/t, the result reflects the evolution of the 
corresponding spin-spin correlator in the Heisenberg model, 
Fig. [TS] Finite-size effects have been corrected and first order 
hopping contributions entering through the nearest neighbor 
correlator {cL\aj) x were removed (see text). 
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Jt 

FIG. 15: Evolution of the momentum space density-density 
correlator X]fe('^T'=+T'^Tfe)^ foi" ~ (minus the trivial 
parts on the right hand side of Eq. (|25p ) in the Hubbard model 
with U = 8 and A = 16, 24. The Hubbard results follow once 
more the Heisenberg curves, except for some quick oscillations 
due to first order hopping processes which die out for |A| far 
from U . 



Integrable many-particle systems do not relax to the 
well known canonical, or "thermal" , ensembles (a fact 
that was already observed experimentally e.g. in [5l[). 
If they relax the steady state is due to the integrals of 
motion to a much more constrained ensemble [2^, Isot . 
This could be detected experimentally by comparing the 
steady state correlation functions after time evolution to 
those obtained for the corresponding thermal ensemble. 
The temperature should be chosen such as to have the 
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FIG. 16: The same observable as in Fig. 1141 except, that the 
observable has not been corrected for finite size effects and 
the effects of the correlator {a\aj)^- We see here clearly, 
that those corrections of the raw data are important to 
achieve comparability to the corresponding Heisenberg result 
in Fig. [la 

same energy in both states. For (free) integrable mod- 
els, the relaxation occurs due to a phase averaging ("de- 
phasing" ) effect [1^ . In Section IV Bl the relaxation in 
the Bethe ansatz integrable Heisenberg model is treated 
within a mean field approximation. Also in this case, 
relaxation is connected to a phase averaging effect. 

Nonintegrable systems are generally believed to re- 
lax to a thermal ensemble due to effective scattering 
processes. Recent numerical analysis of such systems 
m, m, HI] is not yet fully conclusive due to limita- 
tions on maximum observation times (density-matrix 
renormalization-group) or system size (exact diagonal- 
ization). Analytical approaches are usually restricted to 
rather exotic models or limiting cases. See e.g. [13,111] for 
investigations by dynamical mean- field theory (DMFT). 

In our setup, the nonintegrable two-species Bose- 
Hubbard model could be tuned so close to the Heisenberg 
regime (large ] [/^ — ] ) that thermalization occurs very 
slowly. One might hence observe first a relaxation to a 
nonthermal (almost) steady state due to the integrable 
Heisenberg dynamics, which would then be followed by 
slower thermalization due to the remaining nonintegrable 
first order processes of the full Bose-Hubbard Hamilto- 
nian. 



B. Relaxation for the Heisenberg magnet in mean 
field approximation 

In this section, we investigate analytically the relax- 
ation of the Heisenberg magnet with the initial state be- 
ing the Neel state (HH). In particular we will derive that 
the (staggered) magnetization decays as 1/i'^/^ due to a 
phase averaging effect. 

The model is Bethe ansatz integrable [1^ How- 
ever, it is in general not possible to solve the equations of 
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FIG. 17: For the initial state (|14|) . evolution of the entangle- 
ment entropy with respect to a partition of the system into 
left and right half. The growth is roughly linear in time (com- 
pare e.g. to [s^]) resulting in an exponential increase in the 
computation time required for the simulation. The more im- 
portant first order processes are, the faster the entanglement 
entropy increases. The tuning from the Heisenberg model 
(jC/^ — A^] oo), where no first order processes occur, to 
the regime ]A] ~ (7 = 8 can be understood as a smooth 
increase in the number of relevant degrees of freedom, result- 
ing in a stronger entanglement growth. The entanglement 
entropies for the Heisenberg ferro- and antiferromagnet are 
identical, because the corresponding density matrices are in 
the {Sf }-eigenbasis simply related by complex conjugation; 
cf. Section HVBl 

motion for arbitrary initial states. With appreciable nu- 
merical effort this has been achieved recently (only) for 
the initial state being the groundstate plus a one-particle 
excitation [s^ . To investigate the dynamics neverthe- 
less, we hence employ a mean field approximation for the 
Sf+i term. 

^ = {\{StSr^, + 5-5+ i) + 5f 5f+i) 

- {\{StS-+, + 5-5+1) - 2(-l)V.(i)5f ) (27) 

where the order parameter is the staggered magneti- 
zation 

p^^T^ E.(-ir(^^) = ^ E.(-irK " \)- (28) 

After a Fourier and a Jordan- Wigner transformation (ssl . 

m with a := (-l)^"='i(^'.+5)5- and 5f = cjq - i, the 
mean field Hamiltonian and the staggered magnetization 
read with £fc := cos A: 

(£fc4c, - 2p.(t)4+^c,), (29) 

-f <fc<f 

= ^ E 25R(cic,^J. (30) 
The initial state is the Nccl state and reads in the 
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fcrmionic operators for t = with Uk{0) = Wfc(O) = l/\/2 
\m= n iMt)cl+Vkit)cl^J\0). (31) 

— -2<fc<-2 

So each mode Cfc is in the initial state only correlated with 
mode Ck+TT- As the mean field Hamiltonian ([29|) couples 
for every k also just those two modes, the state remains 
in the form ([?T|) for all times. With ihdtCk{t) = [ck, H{t)] 
one obtains the equations of motion (h = 1) 

idtUk{t) = Ek ■ Uk{t) - 2p^{t) ■ Vk{t), (32a) 
idtvk{t) = -Ek ■ Vk{t) - 2p^(t) • Uk{t), (32b) 

a system of N coupled nonlinear differential equations. 

Those can be integrated numerically, yielding for the 
staggered magnetization ^^(i) a damped oscillation de- 
caying as l/t^/'^. Fig. [TH] compares p-K{t) from the mean 
field analysis to the corresponding DMRG result (Fig. [3]) 
and shows good qualitative agreement. 




FIG. 18: Evolution of the magnetization on a particular site 
X, starting from the Neel state and evolving with respect to 
the isotropic Heisenberg Hamiltonian, once with DMRG and 
once in the mean field approximation. The (staggered) mag- 
netization shows a decay (blue). In the mean field ap- 
proach one sees that local relaxation is connected to a phase 
averaging effect as is typical for integrable models [2y]; see 
text. 

As demonstrated in [sl], where the same equations of 
motion were obtained for the evolution of a system of 
spinless fermions, Eq. (|32|) is equivalent to the equations 
of motion of the classical Hamiltonian 

Hs = - '^^f^Sl + ^Y.^SlSi+SlSD (33) 

— 7T<.k<7T k^k' 

with the Anderson pseudospm variables 5*^ = v'^.Uk, 
SI = \{\vk? - iMl'Ifor -f < < f and Sl_,, = S't, 
^k+TT = ""^fc'^i This Hamiltonian occurred in 

the mean field analysis of quenches in fermionic conden- 
sates; see e.g. [13, [H, [s^ [60] . From this, it is known that 



(|33)) and hence JSH) are integrable ^ due to the N/2 
integrals of motion L^, 

ifc^e, + 2V^^, dtLl = 0. (34) 
^k' ^'^ - 

One can now argue that the x and y components of the 
vectors Lk will vanish for large times, as done in [s^ . 
From this one can determine the (nonthermal) steady 
state, by equating {Lf.{t — > oo))^ with L|(t = 0). The 
resuh is limt_oo 5*^ (*) = ^cosA:. With 5"^ = ^(KP - 
|ufcp) and l-Ufcp -I- \uk\'^ = 1 Vj, it follows that 

lim \uk{t)\ = VT~cosk/V2, (35) 

t — >oo 

lim \vkit)\ = Vl + cosfc/V2. (36) 

t — ^oo 

With the knowledge of the steady state, the Xji^l"^ decay 
of the magnetization can now be derived. 

To this purpose let us first recall the general dephas- 
ing scenario for d-dimensional (free) integrable models. 
In it was demonstrated that local observables G(t) 
(i.e. correlators) lead in general to expressions of the form 

G(t) = Go + / d''fce'''^('=)7(fc), (37) 

where the amplitude j(h) is determined by the chosen 
observable, the initial state, and the eigenbasis of the 
Hamiltonian. The phase function (/5(fc) is determined by 
the spectrum of the Hamiltonian. Now, the quantity G(t) 
relaxes to Go for large times if the phase function varies 
quickly enough in regions of the h space where the am- 
plitude /(fc) is nonzero. Whether and how quickly an 
observable relaxes is in particular determined by contri- 
butions from points where V'(fc) is stationary or /(fe) di- 
verges. For the paradigmatic scenario of </j(fc) ~ (^o + |fc|^, 
/(fc) ~ l/fc™ near a stationary point fco = 0, the integral 
in (P7|) behaves as 

givot J d'^fc |^e*l'=l'* - / dg^e*?*, x = ^^df^. (38) 

Hence the time-dependent contribution to G{t), for t 
oo, does not vanish if x ^ Ij vanishes as l/t^~^ if < 
X < 1, and at least as 1/t if x < 0, [75l |. 

Now we come back to the staggered magnetization. 
Expressed in the variables u and w, it reads after going 
to the thermodynamic limit 

p^{t)^ jdm{ulvk)^ y"dfc3?(e*'^WV(fc,i)). (39) 

This is, except for the additional time dependence of 
the amplitude function f{k,t), an integral of the form 
P7p . Presuming that vanishes for long times, it fol- 
lows from the equations of motion ([5^ that for large t, 
the phases of Uk and Vk are roughly zLsk ■ t and hence 
Lp{k) ~ Sk — £k+Tr = 2cosfc, which is stationary (with 
£ = 2) at fc = 0; cf. Fig. [TH For finite times, the 
amplitudes of Uk and Vk are finite as \uk\'^ = (nk) and 
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|ffc|^ = (jT-fc+Tr) = 1 — ("^fe)- Hence m = for fc = which 
is also confirmed numerically in Fig. 1191 The dephasing 
of the staggered magnetization ((39|) is determined by the 
stationary point k = of ip. With d = 1, m = 0, and 
£ = 2 we have x ~ ™+^^~'^ — i. The phase averaging 
accounts hence for a factor of l/i^"'^ ~ l/t^/"^ for the 
decay of the staggered magnetization. Linearizing the 
equations of motion around the steady state we find that 
IuqUqI {f{k,t) around fc = 0) decays as l/\/t and further 
that only a vicinity < l/v^offc = Ois contributing 
to (the leading order of) the integral ((39)l . Taking all 
this together, we infer the l/t'^/^ decay of the staggered 
magnetization. 




FIG. 19: Evolution of the occupation number (uk) = 1 — 
{Uk+n), phase (p{k) = arg{ul{t)vk{t)) - aTg{ul/2i't)'^^/2{t)) of 
the wavefunction, and |u^i!(i| for each pair of modes Ck, Ck+n- 
The initial state, the Neel state with Uk = life — is 
evolved with the mean field approximation of the isotropic 
Heisenberg Hamiltonian (|29|) . Except for the lowest panel, 
the curves for t — 128 coincide (within resolution of the plots) 
with the limiting curves for t oo which are ln^li^ 2cosfc, 
and I sin fc 1/2, respectively. 



VI. VALIDITY OF THE EFFECTIVE SPIN 
MODEL 

One may be wondering why the restriction to the 
single-occupancy space Hi, ([5]), is justified (if the ini- 
tial state of the system is in Hi and we evolve with the 
effective Hamiltonian) , although the coupling to the rest 
of the Hilbert space has the same strength as the coupling 
for dynamics inside the subspace Hi and although parts 
of the rest of the Hilbert space overlap energetically with 
Hi- There could be considerable transition rates out of 



the subspace with (predominantly) one particle per site, 
rendering a description or comparison with dynamics of 
the effective model derived for that subspace useless. We 
will assess here that this is not the case (for the large-fJ 
limit). 

First of all. the numerical results of Section ITVl showed 
that for large — A^|, the Hubbard curves follow 
quite precisely the Heisenberg curves, indicating very lit- 
tle transitions to other subspaces. One can also give a 
somewhat hand waving but rather suggestive argument. 
We will show in the following that transition matrix el- 
ements leading out of Hi occur predominantly to states 
with energy difference ~ and diminish in the large-[/ 
limit. Those yield therefore finite small transition ampli- 
tudes. In higher orders of the perturbation theory, there 
are also transitions to states with energy ~ f7± A, which 
will lead to a small (controllable) transition rate out of 
Hi. 

In the full effective Hamiltonian (jA16p . we regard the 
term V = i[S,H^] that generates or destroys double 
occupancies, i.e. generates transitions between subspace 
Al" with n doubly occupied sites perturbation. 

= + (40) 

The subspaces Al" separate (energetically) further into 
Mm-fj, m with m doubly occupied and fi empty sites on 
sublattice A {n — m doubly occupied and n — fi empty 
sites on sublattice B); i.e. Hi = q. Fig. [^D] shows the 

many-particle spectrum of the effective Hamiltonian H^ff 
for the subspaces A^" and A^^ as obtained from exact 
diagonalization in the 5^ = sector for = 8 sites. 
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eigenenergy E/t 

FIG. 20: The many-particle spectrum of the effective Hamil- 
tonian H^ff for the subspaces Af" (green) and (black) 
with f/ = 8 and A = 2, 4, 6, 8, 10, V^U, 14, 16, . . . as obtained 
from exact diagonalization in the — sector for N = 8 
sites (dimAl" = 70, dimM^ = 2800). Each dot corresponds 
to an eigenenergy. For the plot, small random numbers were 
added to the A values to give a rough impression of the den- 
sity of states. The subspace Al^ is separated energetically 
into Mq U Mo i around E = U and Ml i, M-i o around 
E = U±A. ' 
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The single (quasi-)particle excitations in these sub- 
spaces have energies of order 0{J,t) - spinwaves and 
hopping of doubly occupied and empty sites. However, 
the subspaces overlap energetically (in the thermody- 
namic limit) as, in a qualitative picture, one can have 
^ N quasi-particle excitations resulting in the width ~ 
A^|J| ^ U of the spectrum for each subspacc. Specifically 
for A^°, the lower and upper bounds on the spectrum 
are determined by the ground state energies of the ferro- 
magnetic and the antiferromagnetic Heisenberg models. 
Those are in the thermodynamic limit Ef^ = —jJN and 
-£;ato = (ln2~i)J7V, [3i,[6li. 

If we act on a state \tp) S of energy E with the 
operator V = i[S,H^] (cf. Appendix [X)) . firstly, S gener- 
ates a double occupancy and an empty site; |(Ti, (Ti+i) i-^ 
|(TiiTi+i , 0) on two neighboring sites i and i -|- 1. Secondly, 
a corresponding hopping term from Hi acts on i (or i + 
and i — 1 (or i + 2) such that e.g. |(Ti_i, difTi-i-i, 0) i-^ 

|o-i_lCrj, (Tj+l, 0) or |(Ti(Ti+l, 0, 0-4+2) IfifTj+l, (Ti_|_2, 0). 

Hence both, the doubly occupied and the empty site are 
on the same sublattice ^ or ,S and the resulting state 
V\tp) gMIqU Mil has energy E ^ E + U: see Fig. HH 
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FIG. 21: The subspace , with exactly one doubly occupied 
and one empty site separates energetically into three different 
subbands Ml^i, Mo,o^Ml^i, and AlLi,o- The operator V = 
i[S,Ht] maps states from Ad" to states from TMoq U A^o.i 
that differ in energy by ~ f/ (see text). 

Let us consider transitions from to M^. For any 
initial cigcnstate \i) € the transition amplitude, to a 
state I/) € Al^ is in the Born approximation given by 



destroy in a qualitative picture only a small number of 
quasi-particle excitations as it is a product of only four 
ladder operators. This will change the energy only by a 
small amount of 0{J,t). So hcufi will indeed be of order 
0{U) for all nonvanishing transition amplitudes {f\V\i). 

To illustrate this. Figs. and E51 show the transition 
matrix elements {f\V\i) between eigenstates of the effec- 
tive Hamiltonian H^^ for the subspaces Al" and Ai^ a.s 
obtained from exact diagonalization in the 5'^ = sector 
for = 8 sites. 
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FIG. 22: Transition matrix elements {f\V\i) between eigen- 
states of the effective Hamiltonian H^g for the subspaces Al" 
and A4^ with [7 = 8 and A = 16 as obtained from exact 
diagonalization in the — sector for N — 8 sites. Each 
dot corresponds to a nonzero transition matrix element. The 
narrow panels to the left and bottom show the corresponding 
eigenenergies. Nonvanishing matrix elements exist only for 
states with energy difference of 0(^7). 

Small matrix elements to states with energy difference 
[/ ± A remain. For unfortunate choice of U and A one 
may hence encounter nonvanishing transitions to states 
in A^">o with hojf, ~ 0. Consider e.g. A = 2U. In 
this case, two actions of the operator V may end up in a 
state I/) with comparable energy {U + {U — A) = 0) and 
hence to a (finite but small) transition rate out of Ai^. 
By appropriate choice of the ratio U/A, one can achieve 
that the effect occurs only in higher orders V, resulting 
in a small transition rate. Further, the transition matrix 
elements itself can be made small by going to the large-[/ 
limit dZl). 



This estimate of a small (oscillating) transition amplitude 
follows from the consideration that nonvanishing transi- 
tion elements exist only for states with energy differences 
huifi = 0{U). We have pointed out that the subspaces 
and ultimately overlap energetically. However, 
states I/) and \i) with comparable energy will have a 
vanishing transition matrix element: As argued above, 
the operator V generates states from Ai^ and causes a 
change of [/ in energy. Besides this it can create or 



VII. PREPARATION OF THE 
ANTOFERROMAGNETIC GROUNDSTATE BY 
ADIABATIC EVOLUTION 

In the Sections llVl and lVll we have given arguments and 
gathered numerical support for the fact that transition 
rates from the single-occupancy subspace 7i°"^, (fTT|) . to 
the rest of the Hilbcrt space can be made small for time 
evolution with the Hubbard Hamiltonian. If this can 
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FIG. 23: Transition matrix elements (/IV'li) between eigen- 
states of the effective Hamiltonian for the subspaces 
and with L*' = 8 and A = 10 as obtained from exact 
diagonalization in the 5*^ = sector for N = 8 sites. Still 
nonvanishing matrix elements exist only for states with en- 
ergy difference of 0{U). But as A is closer to U here, the 
matrix elements are larger in amplitude, (|41|l . and the spec- 
tral subbands are broader due to a larger effective coupling J. 



also be realized experimentally for sufRciently long times, 
it would be possible to prepare for example the ground 
state of the antiferromagnetic Heisenberg model by adi- 
abatically switching on the coupling t' between initially 
isolated double wells. Fig. [I] i.e. switching from t' = 
to t' — t, while t is kept constant. As demonstrated in 
[25l . [26j for the initial situation of isolated double wells, 
the groundstate of the single occupancy subspace ([5]) can 
be prepared experimentally. 

For the adiabatic approximation [6^ . [63j to be appli- 
cable, the system needs to be gapped on the whole path 
in the space of system parameters except for the end 
point, where the gap has to close abruptly enough. As 
argued in Section IVTl transitions to other subspaces with 
(quasiparticle) double occupancies can be neglected for a 
certain period of time T that can be made very large. So 
we only need to worry about transitions from the 7i°"^ 
groundstates to excited states inside the subspace, i.e. 
we need to derive conditions on the dependence of the 
corresponding energy gap on the hopping t' between ini- 
tially isolated double wells such that t' = t can be reached 
adiabatically in a finite amount of time t < T. 

The quantitative condition for adiabaticity is generally 
stated as 



{Eo{t)\^\E„it)) 



<1 V, 



t£[0,T],n^Oi 



(42) 



where \En{t)) label the energy eigenstates and \Eo{0)) 
is the initial state. Recently, substantial problems were 



pointed out [6J, |65| and two more conditions added [66 



dt 



d {E,{t)\^\EUt)) 



dt Eo{t)-En{t) 



« 1, 



(43) 



dt 



{E,{t)\ 



^\En{t)) 



Eoit)-E„it) 



\{Er.{t)\^\E^{t))\^l. 



(44) 



If we have one time-dependent system parameter p{t), 
namely the dimerization p = S, where 



S :-- 



1 - \J'/J\ 1 - |t'/tp 



i + i + W 



(45) 



and one part of the Hamiltonian is linear in that param- 
eter (this is the case for the effective Hamiltonian and 
S ^ 0), the numerators of are proportional to 

the sweeping speed v{t) :— dp(t)/dt. The denominator 
is the spectral gap Eg{t). Only points pir) in parameter 
space where the gap vanishes are problematic. If the gap 
vanishes as 



^g(t)« b(r)-p(t)r, iy>0, 

V should (for t close to r) be reduced as c|r — 
cording to 



1 > 



Egit) |/;"*ds.sA'|- 



fi-v{fj.+ l) 



(46) 
Ac- 

(47) 



Hence, only for ly < 1, i.e. for gaps that close abruptly 
enough, adiabaticity can be reached with fi > The 
second condition. (l43l). is in this scenario fulfilled auto- 



matically, the third, (jH]), implies ^ > which is also 
true. 

For J' ~ J, a situation which was examined inten- 
sively in the context of spin-Peierls systems, the model 
was first treated by Jordan- Wigner transformation and 
subsequent bosonization [67j . The precise result for the 
excitation gap can be obtained by a mapping to the four- 
state Potts model or conformal field theory [6^ (see 
also [Z^])- The gap is given by 



ln(5|i/2 



(48) 



That means we have a gap with v — 2/3 < 1 and 
hence the gap can be closed in a finite amount of time 
with exponent /i = 2. This means that the dimerization 6 
has to be varied with speed v{t) = c|t — <p, hence S{t) = 
f |t - tp. One needs thus the time t = (2/c)^/3. The 
smaller c is, the farther we arc in the adiabatic regime 
but the longer we need for the preparation. An analysis 
of how small c is to be chosen to achieve a given accuracy 
of the prepared state could be carried out along the lines 
ofRef. [nl- 

Note that in [t^, it was recently discussed within a 
mean-field approach, how the antiferromagnetic phase 
of the three-dimensional Fermi- Hubbard model could be 
reached by adiabatic tuning of the lattice potential. 
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VIII. CONCLUSION 

We have studied a setup of two species of ultracold 
bosonic atoms in an optical superlattice, which realizes 
in a certain parameter regime the Heisenberg ferro- and 
antifcrromagnet. The focus was in particular on time 
evolution of nonequilibrium states. Our numerical results 
and analytical considerations showed that the physics of 
Bose-Hubbard model implemented in the experiment dif- 
fers for certain parameter ranges considerably from the 
physics of the effective Heisenberg models. Note that 
this would also be true for alternative suggestions as in 
0, H^j m, ■ The spin states up and down can in 
general not be identified directly with a bosonic particle 
of one specific species. The regime where the correspon- 
dence between the two models is good, implies higher 
requirements on cooling and coherence (coherence time) 
in an experimental realization. The explicit form of the 
Schrieffer- Wolff transformation was used to analyze the 
transition rates out of the magnetic subspace of the full 
Hilbert space. 

In contrast to the accomplished experiments [1^ [2^ 
for isolated double- wells (filled each with two particles), 
the setup of coupled double-wells discussed here allows 
for relaxation of the many-particle state. In the numerics 
we observed indications for (local) relaxation to steady 
states. For the Heisenberg model in a mean field approxi- 
mation, we explained how the relaxation is connected to a 
phase averaging effect. This is typical for integrable mod- 
els which have nonthermal steady states. Nonintegrable 
models are generally believed to thermalize due to effec- 
tive scattering effects. Our setup can be tuned from the 
nonintegrable Bose-Hubbard model to the Bethe ansatz 
integrable Heisenberg model and could hence be used to 
study the differences of the relaxation processes experi- 
mentally. 

Finally we argued that the groundstate of the Heisen- 
berg antifcrromagnet could be prepared by tuning an al- 
ternating hopping parameter of the superlattice adiabat- 
ically. 
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APPENDIX A: DERIVATION OF THE 
EFFECTIVE MODEL BY SCHRIEFFER- WOLFF 
TRANSFORMATION 



Here we derive the effective spin Hamiltonian de- 
scribing the physics of the two species Bose-Hubbard 
model dll) in the subspace 7i°"^, (fTT|) . where every site is 



occupied by exactly one quasi-particle. The spin-spin in- 
teraction is generated by second order hopping processes 
of the particles. In the large-C/ limit, transitions from 
H™^^ to bands with double-occupancies are energetically 
hindered. In particular, transitions from 

Til = span{|ni, . . . ,nAr),n|i + nn = 1 VJ (Al) 

to the subspace with double-occupancies (and holes), 7^2, 
can be treated perturbatively. 

The Hamiltonian contains terms, linear in the hop- 
ping t, which couple Tii to the subspace with double- 
occupancies. We are looking for a canonical transfor- 
mation H ^off" ■— e.'^^He~^^ such that the single- 
occupancy subspace TLi of the resulting quasi-particles 
e^^a^ie^^^ couples only in second order to double occu- 
pancies (7^2)- 

The calculation can be done in analogy to the deriva- 
tion of the Kondo lattice model [t^ from the periodic 
Anderson model [H, |4l|, [tI] or the t-J model [lH from 
the fcrmionic Hubbard model and is modified only 
by the asymmetry term and the finite intra-species 
repulsion (double occupancies | tt) and | [[))■ 

Let us rewrite the Hamiltonian ([1]), restricted to the 
subspace Tii U 7^2; in the form 



Hq = Ha + Hjj 

Ha = ^ '^iTT-o-A 



Hu = U'^n^iUii + ^ ^ n„i{n„i - 1) 



A± = -t^(f4+t%), 



(A2) 
(A3) 
(A4) 

(A5) 
(A6) 
(A7) 



where (ij) runs over nearest neighbors (one index from 
each sublattice), H^ increases/decreases the number 
of doubly occupied sites and i/^ leaves it unchanged; 
Fig. [21 

(A8) 
(A9) 
(AlO) 
(All) 

where S denotes the Kronecker delta and in its argument, 
rii = niji + nil denote the particle number operators. 
For the operators T^^ , we further distinguish between 
those which change the number of (a) inter-species and 
(b) intra-species double occupancies, see Fig. [Ml 
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FIG. 24; Illustration of the hopping terms (|M)) - (|MT|) . T*-^ 
increases/decreases the number of doubly occupied sites and 
T"''^ leaves it unchanged. 



1. Schrieffer- Wolff transformation 

The Haniiltonian contains terms , linear in the hop- 
ping t, which couple Hi to the double-occupancy sub- 
space. We are looking for a canonical transformation 
H — > H^g" = e^^ He'"^^ such that the single-occupancy 
subspace Hi of the resulting quasi-particles e^^a„ie~^^ 
couples only in second order to double occupancies (7^2)- 

= H + i[S,H] + - [S, [5, H]] + 0{S^H) 

= Ho + + H+ + + i[S, H] + -[5, [S, H]] + . 



In the first commutator, the contribution [S, Hq\ domi- 
nates and we therefore look for a generator S such that 



(A13) 



From this equation follows with iJo = 0{U,Ust^) and 



Hi = 0(1) that S = 0{ 



) and hence 



2. Effective spin Hamiltonian for half filling 

The full effective Hamiltonian reads 

-^cff ' = ^i + Hu + Ha 

+ '-[S, H+ + H^] + i[S, H^] + 0(f3). (A16) 

The commutator terms still couple Hi with the rest of 
the Hilbert space (subspaces with differing numbers of 
doubly occupied sites). However, this coupling is now not 
0(f), as in the original Hamiltonian jAS]), but of O(t^). 
This was achieved by the Schrieffer- Wolff transformation 
S, which replaces our original particles aa-i, by particles 
with a cloud of hole-double-occupancy fluctuations a^i — > 
e^^ aa-ie~^^ . In the single-occupancy subspace Hi at half 
filling, (|XT|), H^, Ha, Hu, i[S,H^] and the terms of 
third order in the hopping are all ineffective such that we 
arc left with 

iJeff HX\n, = ^[S, Ht + + (A17) 

The commutator consists of hopping terms via virtual 



-7 

and can be rephrased as spin-spin interactions. With 



double-occupancy states. They are of the form T^i-^T^^^ 



r-,?t%'k. = 1 + ^StS^, (A19) 

Erp~a rp+a\ Q+ Q— I Q~ Q+ 

the effective Hamiltonian (jAl7P reads 



(A18) 
(A19) 
(A20) 



H, 



cff - 



(ij) 



■ sfsy) 



on'), (A21) 



where (see also [2g|) 



HX = Ho+H^+i[S,H,] + -[S, [S,H]]+Oi 



(f/,C/.,A)2' 



This yields for (|A13[) the solution 
t 



rp+a 



+a 



i 



-T, 



+b 



t 



Us + Aj - Ai 



h.c. 



(A14) 



With Ai from we can now state more precisely 

t t 



s = o 



U ±A Us±A 



(A15) 



So the perturbative treatment will break down near the 
crossing point A = U and for large hopping t. 



41-^ J7 

[/2 - A2 



Js = 2- 



(A22) 



This Hamiltonian is, except for higher order effects, the 
XXZ model. In the bulk of the article we specialize to 
U = Us, i.e. J — Js = — J and have hence the isotropic 
Heisenberg ferromagnet for A < J7 (J > 0) and the 
isotropic antifcrromagnct for A > [/ (J < 0). As was 
already pointed out, the full effective Hamiltonian (|A16[) 
still contains a coupling to the subspace with one double- 
occupancy (of quasi-particles). The approximation made 
by neglecting it is discussed in Section fVll A peculiarity 
of our situation is that, due to half- filling of both parti- 
cle species, we are restricted to the = sector of the 
Heisenberg model. 

One can go to higher orders in the perturbative treat- 
ment of the hopping, by adding higher order terms to the 
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generator S of the Schrieffer- Wolff transformation. The 
next order term 5^^^ has to be chosen such that it eUmi- 
nates the term i[S'^^\ H^] in (jA17p . This would result in 
a further contribution to the effective spin Hamiltonian, 
namely next nearest neighbor and four-spin interactions, 
generated by sequences of four virtual hopping events. 

APPENDIX B: POSTPROCESSING OF 
DENSITY-DENSITY CORRELATORS 

1. Elimination of finite-size effects for the numerics 

The numerics were done for a finite size system (cf. Sec- 
tion |TVC|. To correct for the resulting finite-size effect is 
simple for the calculation of the spin-spin correlators in 
the Heisenberg model in the r.h.s. of (|23p and p4)) . One 
can use 

(S'f S'j^)^ H-> (5'"^5'f^_i_^._,.)^, (Bl) 

where Xi is some site in the middle of the system that 
is odd (even) for odd (even) i. This corresponds to the 
invariance of the systems under translations by multiples 
of two sites in the thermodynamic limit. 

The momentum-space density-density correlators, 
(|25p , in the Hubbard model are determined from the real- 
space four point correlators ([23ll . To eliminate finite-size 
effects of those, we note first that due to the restriction 
of all correlations to a (causal) time-space cone, the four 
point correlators behave for large distances as (spin in- 
dices suppressed) 

Cijn {O'i'^jO'n+j-iO'n) 

+ [(«!««) + iS^.n - S^,j){n,)]. (B2) 
Consequently, the quantity 

- [(«I«n)(«i+,-««j) + i^n - <5,,,)(n,)] (B3) 

is localized; it has support only for j and n inside the 
causal cone centered at site i (cmp. to Section lIVEj) . So 



the value of Cijn for the thermodynamic limit is approx- 
imated well by 

^ijn '■~ ^ijn "t" 9i,j9n+j — i,n 

+ [9i,n9n+j-i,j + iSi,n - Sij)gis], (B4) 

where gi j is the (approximate) single particle Green's 
function in the thermodynamic limit 

9i,j = {al,a^,+j-t), (B5) 

and for odd (even) i, Xi is an odd (even) site in the 
middle of the system. 

2. Reduction of effects from first order processes 
for numerics and experiments 

As disclosed by equation (|B2p or (|B4[) the four point 
correlators entering the momentum-space density-density 
correlator contain contributions from products of single- 
particle correlators {a\a^) ^. Those are trivial when evolv- 
ing with the Heisenberg model: {a\aj)^ = 6ijn^i{t), as 
there is exactly one particle per site. But according to 
(IT71) . they have contributions of 0(5^), when evolving 
with the Hubbard Hamiltonian. In the comparison of 
observables evolved with both models, those correlators 
enter hence as a major carrier of disturbance. To achieve 
comparability it would be desirable to remove contribu- 
tions from (aja^) completely. This would be possible for 
our numerical analysis. In a corresponding experiment 
however, the quantities are not available. Hence we con- 
fined ourselves to removing only the contributions from 
nearest neighbor correlators (ala^.^.^). That means gi,j in 
(|B4p is set to zero for j = izL 1. This was already suffi- 
cient to demonstrate the correspondence of the dynamics 
if we are safely in the large-t/ limit ([7]). An experimental 
procedure for the measurement of the nearest-neighbor 
correlator was suggested in [4^. Hence, the same manip- 
ulations might be carried out for experimentally obtained 
momentum-space density-density correlators. 
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